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ABSTRACT 

An  investigation  into  the  problems  of  economically  pro- 
cessing sewage  on  board  Naval  ships  resulted  in  the  develop- 
ment of  two  computer  simulations  employing  Monte  Carlo 
analysis  to  describe  the  generation  of  sewage.   Simulation 
one  was  based  on  a  non-homogeneous  Poisson  process.   For 
simulation  two,  an  empirical  distribution  describing  the 
arrival  behavior  of  sewage  to  the  holding/processing  unit 
over  a  24  hour  period  was  applied  to  known  data  on  sewage 
generation . 

Results  of  the  two  simulations  were  compatible  with  one 
another.   Aside  from  pointing  out  a  most  feasible  combina- 
tion of  holding  tank  capacity,  processor  rate  and  processing 
policy,  the  simulations  also  indicated  that  a  revision  of 
the  Navy's  design  parameter  for  the  daily  per  capita  sewage 
generation  rate  was  in  order.   The  simulations  were  designed 
so  that  use  by  others  with  their  own  data  would  be  an  easy 
matter. 
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I.   INTRODUCTION 

A .   BACKGROUND 

During  the  time  between  the  conception  of  this  thesis 
and  its  completion,  significant  changes  in  both  policy  and 
attitude  have  been  made  by  the  Navy  in  the  field  of  pollu- 
tion abatement.   These  alterations  have  resulted  from  in- 
creasing public  pressure  for  the  ecological  restoration  of 
U.S.  inland  waters  and  subsequent  congressional  legislation 
which  outlined  future  restrictions  to  all  seagoing  vessels 
and  watercraft.   More  specifically,  in  August  19  71,  the 
Chief  of  Naval  Operations  stated  "that  the  water  pollution 
abatement  program  for  Naval  ships  has  as  its  long  term  goal 
the  installation  of  two  basic  systems  aboard  ship: 

1.  A  Marine  Sanitation  Device  (MSD)  which  processes 
sewage  to  an  acceptable  effluent; 

2.  A  Collection,  Holding  and  Transfer  (CHT)  System 
which  is  designed  to  transfer  all  liquid  wastes  (except 
oil)  ashore  or  to  a  barge." 

Late  in  1971,  however,  the  Senate  issued  further  re- 
strictions to  vessels  to  the  extent  that,  by  the  year  19  81, 
there  shall  be  no  toxic  effluent  discharge,  whatsoever, 
from  any  ship  into  U.S.  inland  waters.   As  a  result,  Naval 
design  efforts  refocused  primarily  on  the  development  of 
a  Collection,  Holding  and  Transfer  System.   This  shift  in 
emphasis  was  also  aided  by  the  fact  that  prototypes  of 
Marine  Sanitation  Devices  had  suffered  from  low  system 


reliability  and  excessive  costs.   The  present  goal,  then, 
of  the  Navy  is  the  development  of  a  sewage  holding  tank 
system  with  the  capacity  to  retain  12  hours  of  sewage  gen- 
erated by  a  full  complement  crew  and  with  the  capability  of 
transferring  such  sewage  to  an  external  processing  facility 
ashore.   Such  a  system  is  to  be  installed  in  all  ships  re- 
maining in  the  fleet  after  19  77. 

B.   STATEMENT  OF  THE  PROBLEM 

It  is  the  contention  of  the  author  that  neither  strictly 
holding  sewage  for  later  transfer,  nor  immediately  process- 
ing the  sewage  internally  as  it  arrives  to  the  processing 
unit,  is  the  most  feasible  economic  solution  to  pollution 
abatement.   A  policy  of  holding  sewage  and  transferring  it 
to  some  other  processing  station  would  greatly  hamper  the 
Navy's  capability  to  act  as  a  mobile  self-sustained  unit; 
likewise,  the  policy  of  rapid  processing  aboard  would  re- 
quire expensive  units  of  high  output  and  reliability.   It 
is  believed  that  viable  trade-offs  can  be  realized  through 
a  combination  of  the  two,  employing  a  smaller,  cheaper  pro- 
cessor which  could  operate  at  a  lower  processing  rate,  and 
a  smaller  holding  tank  that  could  absorb  sewage  flow  during 
peak  periods  of  the  day. 

A  second  contention,  that  is  perhaps  more  germaine  to 
the  current  Navy  policy  of  a  twelve  hour  holding  system  with 
subsequent  transfer,  is  the  belief  that  the  presently  ac- 
cepted daily  per  captia  rate  of  30  gal/man/day  of  sewage 


generation  is  over  estimated  [Ref.  5].   The  parameter  was 
determined  from  data  gathered  by  the  Marine  Engineering 
Laboratory  on  the  daily  sewage  generation  rates  per  man  on 
board  four  ships.   By  reasoning  that  "any  treatment  unit 
used  aboard  a  Naval  vessel  must  be  designed  to  operate  with 
a  full  complement  24  hours  per  day,"  the  MEL  Report  [Ref. 
18]  recommended  that  the  rates  measured  be  doubled  for  de- 
sign purposes.   This  resulted  in  a  figure  of  30  gal/man/day. 
Since  the  original  figures  were  obtained  on  a  per  capita 
basis,  and  there  was  little  variation,  it  is  felt  by  the 
author  that  more  accurate  results  could  be  obtained  by 
applying  the  unmodified  data  to  a  model  which  also  considers 
the  traffic  pattern  of  personnel  using  the  facilities  over 
a  period  of  time . 

It  was  the  purpose  of  this  study,  then,  to  investigate 
the  problem  of  sewage  generation  and  processing  aboard  a 
Naval  ship,  and  to  develop  a  model  which  would  best  de- 
scribe such  a  process.   This  model  would  be  useful  in  aiding 
decisions  on  the  selection  of  the  most  feasible  treatment 
system. 

The  approach  taken,  as  presented  in  this  thesis,  was  to 
first  formulate  the  problem  analytically  and  to  compare  it 
with  other  research  performed  in  related  cases.   A  computer 
simulation  was  then  developed,  motivated  by  a  similar  study 
done  for  the  air  lines  regarding  sewage  handling  on  board 
large  jumbo  jets  [Ref.  9].   The  simulation,  its  methodoligies 
aod  how  it  can  be.  used  are  described  in  Chapter  III.   The 
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simulation  provides  a  sensitive  model  which  illustrates  the 
trade-offs  obtainable  through  any  treatment  policy  tested. 
It  also  reinforces  the  contention  that  a  revision  of  sewage 
generation  rates,  combined  with  closer  analysis,  can  lead 
to  savings  in  holding  tank  size. 

It  should  be  noted  that  a  simulation  is  no  stronger  than 
the  basic  assumptions  made  regarding  its  underlying  struc- 
ture, and  that  such  assumptions  should  be  examined  criti- 
cally.  With  this  in  mind,  the  simulations  were  developed 
using  the  best  information  available;  however,  they  were 
also  made  flexible  enough  so  that,  should  a  designer  have 
better  arrival  data  or  other  policies  to  test,  such  informa- 
tion could  be  easily  read  into  the  program  to  give  useful, 
real  world  results . 


II .   PROBLEM  FORMULATION 

A.   THE  BASIC  MODEL 

The  overall  process  of  generating  sewage  over  a  specific 
length  of  time,  storing  it  and/or  processing  it  can  be  de- 
scribed analytically.   Further,  it  can  be  related  to  a 
group  of  similar  problems  that  are  generally  classified  as 
storage  problems ,  examples  of  which  are  the  collection  and 
storage  of  water  in  a  reservoir  subject  to  some  release 
rule,  and  inventory  problems  in  which  demand  is  random. 
Any  such  process  in  which  either  the  input  or  output  is 
random  comes  under  the  category  of  a  stochastic  process. 

Such  is  the  case  with  the  sewage  storage  and  processing 
problem,  in  which  the  input,  X(t) ,  (see  Figure  1)  into  a 
storage  tank  is  random  in  nature.   Here  X(t)  represents  the 
amount  of  sewage  generated  by  the  ship's  company  at  discrete 
points  in  time,  t,  corresponding  to,  say,  hours  or  minutes. 
Two  important  characteristics  of  this  input  are  1)  its  time 


X(t) 


Input 


Tank 


Processor 


Figure  1.   Schematic  of  System 
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dependent  arrival  behavior  to  the  tank,  and  2)  its  overall 
distribution.   It  is  assumed  that  X(t)  is  independent  and 
identically  distributed  for  different  values  of  t.   For  this 
discussion,  X(t)  is  assumed  to  be  solely  bodily  waste;  how- 
ever the  model  could  also  be  expanded  to  include  "hotel" 
and  galley  waste. 

Z(t).  represents  the  total  amount  of  sewage  in  the  stor- 
age tank  also  at  time  t.   Aside  from  being  a  function  of  the 
input  random  variable,  X(t) ,  it  is  also  a  function  of  the 
rate,  a,  at  which  the  sewage  is  processed,  and  the  particu- 
lar "release  rule"  or  processing  policy,  ip ,    that  is  applied. 
The  processing  policy  is  essentially  a  switching  function 
which  turns  the  processor  on  or  off  based  on  either  a  time 
plan  (e.g.  on  at  certain  hours  of  the  day) ,  a  tank  level 
plan  (e.g.  turn  on  when  the  level  of  sewage  reaches  a  cer- 
tain height) ,  or  a  combination  of  the  two. 

At  certain  times  of  the  day,  the  processing  rate  may 
exceed  the  input,  X(t) ,  while  at  other  times  it  may  lag  it. 
If  the  processing  rate  is  much  lower  than  the  input  rate, 
Z(t)  will  behave  like  a  sum  of  independent  random  variables. 
For  such  a  case,  the  tank  would  be  expected  to  overflow. 
If  the  processing  rate  is  much  higher  than  the  input,  then 
the  tank  may  never  fill  at  all,  and  Z(t)  will  behave  more 
deterministically  (or  equal  zero) .   One  of  the  questions  to 
be  answered,  then,  is  what  is  the  optimal  processing  rate? 
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In  general,  Z (t)  may  be  described  as  follows: 

K  if  Z(t-1)+X(t) -ctt>K 

Z(t-1)+X(t) -at,  if  0<Z(t-l)+X(t) -at<K 
Z(t)  =  4  (1) 

Z(t-1)+X(t) ,     if  a=0 

,.  0  if  Z(t-1)+X(t)  -at<0 

where  K  is  the  capacity  of  the  tank. 

This  formulation  includes  two  special  cases .   The  first 

is  the  event  that  overflow  occurs  in  which  case : 

r   Z (t-l)+X(t) -at,  if  Z(t-l)+X(t)-at<K 
Z(t)  =  \  (2) 

.  K  if  Z(t-1)+X(t) -at>K 

Here  the  overflow  is  Z (t-1) +x (t) -at-K .   The  second  case  is 

the  event  that  processing  empties  the  tank  completely,  so 

that : 


Z(t-l)+X(t)-at,  if  Z(t-l)+X(t)-at>0 

Z(t)  =   i  (3) 

0  if  Z(t-l)+X(t)-at<0 


B.   OTHER  RELATED  STUDIES 

Equation  (1)  through  (3)  are  not  unique  to  the  sewage 
generation  and  collection  problem.   Rather,  they  parallel 
those  proposed  by  Moran,  Gani ,  and  others  [Refs.  15,  8,  10] 
under  the  broader  discipline  of  the  stochastic  theory  of 
reservoirs  or  storage ,  an  approach  which  has  attracted  con- 
siderable theoretical  interest,  but  has  seen  less  application 

Initially  the  theory,  as  proposed  by'  Moran,  required  the 
inflow,  X(t) ,  to  have  a  known  continuous  probability  dis- 
tribution and  to  be  mutually  independent  of  the  inflow  for 
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any  other  value  of  t  (all  of  which  have  the  same  distribu- 
tion) .   The  level  in  the  tank,  Z(t),  is  subject  to  some  re- 
lease rule,  ty ,  and  the  sequence  (Z(t+1),  t>0}  is  a  Markov 
process,  with  state  transition  probabilities  that  can  be 
computed  from  the  inflow  probability  distribution.   Consid- 
erable effort  has  been  expended  in  finding  solutions  to  the 
steady  state  distribution  of  Z (t)  for  both  continuous  and 
discrete  cases. 

These  initial  studies  were  extended  by  a  number  of 
authors.   One  such  work  by  F.  M.  Lloyd  [Ref.  12]  showed  that 
the  theory  need  not  be  limited  to  inflows  that  were  mutually 
independent,  but  that  it  could  also  be  applied  to  serially 
correlated  inputs  with  the  results  forming  a  bivariate 
Markov  Chain.   Ali  Khan  and  Gani  [Ref.  1]  carried  this  study 
one  step  further  to  consider  the  time  dependent  solution 
to  the  same  model. 

An  excellent  example  of  the  application  of  the  above 
theory  which  is  closely  related  to  the  sewage  problem  is 
the  work  by  Avi-Itzhak  and  Ben-Tuvia  [Ref.  3] .   Here  the 
authors  studied  the  problem  of  "how  big  a  reservoir  to  build 
and  what  should  the  pumping  capacity  be,  such  that  the  cost 
per  cubic  meter  of  water  is  a  minimum."   In  order  to  obtain 
their  solution,  however,  the  authors  first  had  to  assume  an 
approximate  distribution  for  the  daily  input  to  the  reser- 
voir as  well  as  know  the  costs  involved  in  operating  the 
system. 
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In  all  of  the  above  cases  the  underlying  requirement  is 
that  the  distribution  of  X(t)  is  known  or  can  be  approxi- 
mated by  a  sequence  of  discrete  distributions. 
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Ill .   MODELING  THROUGH  SIMULATION 

A.   JUSTIFICATION  OF  SIMULATION 

The  storage  theory  approach  described  above  leads  to 
perhaps  the  most  precise  solution  to  the  sewage  storage 
problem.   It  has,  however,  certain  drawbacks  which  preclude 
its  immediate  implementation.   The  first  of  these  is  the 
lack  of  real  world  information  and  data  with  which  to 
derive  the  necessary  matrix  of  transition  probabilities. 
Secondly,  the  mathematics  involved  are  difficult  and  cum- 
bersome.  This  is  especially  so  when  one  is  attempting  to 
solve  the  problem  for  different  processing  policies  ("re- 
lease rules")  and  processing  rates,  each  of  which  would 
require  its  own  transition  matrix.   Thirdly,  the  complexity 
of  the  distribution  of  X(t) ,  the  input  variable,  whose  sta- 
tionarity  over  time  is  questionable,  would  require  parti- 
tioning the  process  in  order  to  arrive  at  a  solution.   Such 
a  non-homogeneous  process  would  only  further  complicate 
the  analysis.   Lastly,  the  problem  of  a  reversing  boundary 
brought  about  by  abrupt  changes  in  processing  rate,  accord- 
ing to  the  processing  policy  used,  would  add  even  another 
dimension  to  an  already  difficult  problem. 

Taking  the  above  into  consideration,  it  was  decided 
that  a  computer  simulation  involving  Monte  Carlo  analysis 
was  at  least  a  justifiable,  if  not  preferable,  alternative. 
Most  important  to  the  justification  of  a  simulation  was  the 
fact  that  it  could  be  made  amenable  to  any  alterations  in 
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processing  policies  and  processing  rates.   It  is  also  adapt- 
able to  adjustments  in  the  input  distribution  if  and  when 
additional  data  on  the  arrival  patterns  of  sewage  inputs 
are  obtained.   These  characteristics  have  been  incorporated 
into  the  simulations  developed  in  this  study. 

B.   DESCRIPTION  OF  THE  SIMULATIONS  DEVELOPED 

Two  computer  simulations  were  developed  using  FORTRAN  IV 
for  runs  on  the  IBM  360/67  computer.   The  simulations  dif- 
fered only  in  the  method  of  Monte  Carlo  analysis  used  to 
generate  the  sewage  input  distribution,  X(t),  which  was  then 
applied  to  the  basic  mathematical  model  described  in  Chapter 
II  (equations  (1)  -  (3)) .   For  simulation  one,  X(t)  was  de- 
termined by  assuming  a  non-homogeneous  Poisson  process.   For 
simulation  two,  an  empirical  distribution  describing  the 
arrival  behavior  of  sewage  to  the  tank  over  a  24  hour  period 
was  applied  to  known  data  on  sewage  generation.   The  tech- 
niques and  mathematics  involved  are  described  in  subsections 
2  and  3  following  an  overall  description  of  the  basic  simu- 
lation (subsection  1) . 

1.   Basic  Design  Common  to  Both  Simulations 

Aside  from  the  methods  used  to  generate  observations 
of  the  sewage  input,  the  two  simulations  were  identical  in 
the  way  the  sewage  was  processed  from  the  holding  tank. 
(See  general  flow  chart  of  the  simulations,  Figure  2.) 
First  it  was  assumed  that  the  holding  tank  used  in  the 
models  had  sufficient  capacity  so  that  overflow  was  not  a 
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Figure  2.   General  Flow  Chart  of  the  Simulations. 
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problem.   This  was  done  in  order  to  get  some  idea  of  what 
the  maximum  sewage  level  (ZMAX)  in  the  tank  would  be,  ac- 
cording to  the  processing  rate  and  policy  used.   This 
further  simplified  the  basic  model,  so  that  the  only  equa- 
tions needed  for  the  simulation  were: 

fZ(t-l)+X(t)-at,  a>0 

Z(t)  =  {  (4) 

CO  if  Z(t-l)+X(t)-at<0 

Both  simulations  involved  discrete  time  steps;  for 
model  one,  the  unit  of  time  was  5  minute  increments;  for 
model  two,  hourly  increments.   In  both  cases  the  total  period 
of  time  simulated  was  one  week,  based  on  the  assumption  that 
it  was  desirable  to  have  any  sewage  build-up  during  the  week 
processed,  and  the  tank  emptied  by  the  end  of  the  weekend. 
A  new  cycle  would  then  correspond  to  the  start  of  a  new  work 
week.   It  should  be  noted  that  each  simulation  day  started 
at  0500  real  time.   At  each  discrete  time  increment,  a  ran- 
dom observation  of  sewage  input,  X(t) ,  was  generated  and 
the  level,  Z(t) ,  was  computed  using  equation  (4)  for  a 
specific  processing  rate,  a. 

Anyone  of  five  processing  policies,  in  the  form  of 
separate  subroutines  could  be  called  independently  and  ap- 
plied to  the  model.   The  following  five  policies  were  used: 

1)  Plan  1.   Continuous  holding  with  no  processing 
whatsoever  (a=0)  . 

2)  Plan  2.   Continuous  processing  at  some  rate,  a, 
to  be  read  in  by  the  programmer. 
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3)  Plan  3.   Process  continuously  from  0500  until 
1600  each  day;  after  1600  stop  processing  when 
the  tank  level  reaches  50  gal;  continue  pro- 
cessing at  0500  the  next  day.   For  weekends 
stop  processing  anytime  the  tank  level  reaches 
50  gal;  continue  processing  when  the  tank 
level  builds  up  to  500  gal.   (This  is  a  com- 
bination of  a  level  policy.  Plan  4,  and  a  time 
policy,  Plan  5,  as  described  below. 

4)  Plan  4.   (A  level  control  policy)   Stop  pro- 
cessing whenever  the  tank  level  goes  below 
some  prescribed  limit  (50  gallons  was  used  in 
the  simulations  -  the  value  is  read  in  by  the 
programmer) .   Continue  processing  when  the 
level  exceeds  an  upper  limit  (1500  gallons 
was  used) . 

5)  Plan  5.   (A  time  control  policy)   Process  con- 
tinuously for  a  fixed  portion  of  each  day. 
(For  the  simulations  processing  stopped  at 

hour  16,  2100  real  time  and  continued  at  the 
beginning  of  the  next  simulated  day,  0500  real 
time . ) 

Policies  one  and  two  illustrate  the  two  extremes 
possible.   The  policy  of  continuous  holding  (Plan  1)  was  used 
to  determine  the  amount  of  sewage  that  would  build  up  over  a 
certain  period  of  time.   (The  special  case  of  a  twelve  hour 
holding  tank  was  tested  using  this  policy.)   The  policy  of 
continuous  processing  (Plan  2)  was  used  to  get  an  indication 
of  the  minimum  tank  levels  that  could  be  expected  over  any 
range  of  processing  rates  desired. 

The  remaining  three  policies  were  developed  as  ex- 
amples of  probable  real  world  policies.   They  also  served  to 
reflect  the  response  of  the  simulations  to  changes  in  policy, 
for  which  the  simulations  produced  acceptable  results.   Of 
the  three,  policy  five  would  be  expected  to  be  the  least  ex- 
pensive since  it  simply  involves  turning  the  processor  off 


19 


once  a  day.   Policy  four,  on  the  other  hand,  would  require  a 
level  sensing  device  and  some  sort  of  automatic  start/stop 
mechanism  to  the  processor  and  would,  therefore,  entail  a 
considerably  larger  cost.   Policy  three  would  only  require 
manual  operation  and  its  relative  cost  would  fall  somewhere 
in  between  that  of  policies  four  and  five.   The  point  to  be 
considered  for  all  three  is  the  fact  that  the  smaller  the 
time  required  for  the  processor  to  operate,  the  lower  the 
overall  running  cost  and  the  time  available  for  system 
maintenance . 

It  is  emphasized  that  the  simulations  are  in  no  way 
limited  to  just  these  policies  and  that  any  other  policy  may 
be  tested  simply  by  the  addition  of  a  subroutine  similar  to 
those  used  and  a  statement  by  which  to  call  it.   Further 
details  on  the  mechanics  of  using  the  computer  simulations 
are  contained  in  Appendix  D. 

2 .   Input  Distribution  by  Simulation  One 

One  of  the  most  important  problems  to  be  solved  was 
how  to  generate  reasonably  accurate  observation  of  the 
arrivals  of  sewage  to  the  holding  tank.   The  two  parameters 
that  had  to  be  considered  were  1)  the  arrival  rate  and  2) 
the  quantity  of  sewage  per  arrival.   Although  a  study  was 
conducted  in  19  64  [Ref.  19]  to  determine  the  quantity  of 
sewage  generated  per  man  per  day  on  a  Naval  vessel,  no 
definitive  statistics  were  obtained  to  describe  the  arrival 
process  over  the  course  of  a  24  hour  period.   Owing  to  the 
rrature  of  shipboard  routine,  which  can  be  considered 
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reasonably  standard  throughout  the  Navy,  it  was  evident 
that  the  arrival  distribution  would  be  multimodal  over  the 
period  of  the  day.   With  such  changes  in  arrival  rates  it 
was  also  clear  that  the  process  could  not  be  considered 
stationary  for  a  24  hour  period  and  that  some  form  of  non- 
homogeneous  process  would  be  required. 

In  a  related  study  by  Gephart  and  Balachandran 
[Ref.  9]  into  the  reliability  considerations  of  an  airborne 
sewage  system  for  large  jet  liners,  data  was  collected  on 
the  arrival  behavior  of  passengers  to  the  restroom  facili- 
ties.  It  was  concluded  by  the  authors  that  the  arrival 
process  could  be  considered  Poisson  over  "non-overlapping 
and  consecutive  time  segments."   It  appeared  reasonable, 
then,  that  a  non-homogeneous  Poisson  process  could  also  be 
applied  to  the  sewage  generation  problem  on  board  a  Naval 
ship . 

In  its  most  general  form,  the  process  could  also 
be  considered  a  compound  Poisson  process  represented  by 

N(t) 
X(t)  =   Z    Y±  (5) 

i=l 

where  {N(t),  t>0}  is  Poisson  process  and  {Yn ,  n=l,2...}  is 

a  family  of  independent,  identically  distributed  random 

variables  representing  the  quantity  of  sewage  per  arrival. 

If  a  distribution  for  the  body  waste  generated  per  capita 

per  arrival  were  known,  this  could  easily  be  incorporated 

into  the  simulation.   For  the  purpose  of  simplicity, 
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however,  the  quantity  of  sewage  per  arrival  was  taken  to  be 
a  constant  and  the  compound  property  was  disregarded.   Based 
on  studies  [Ref.  19]  which  determined  the  quantity  of  water 
per  flush  to  be  4.5  gallons,  a  rather  conservative  estimate 
of  4.9  gallons  per  arrival  was  used. 

The  final  problem,  then,  was  to  determine  the  best 
estimate  of  the  arrival  behavior  at  head  facilities  over  a 
24  hour  period.   It  was  decided  that  the  day  could  logically 
be  partitioned  into  three  periods :   the  first  being  the 
period  between  0500  and  1600  when  the  majority  of  the  crew 
is  on  board  for  the  working  day;  the  second  being  the  period 
between  "liberty  call"  and  "taps"  (1600-2200)  and  the  last 
including  the  time  from  2200  until  0500.   The  intensity 
function  X(t)  which  describes  the  arrival  rates  could  then 
be  written  for  weekday  periods  as : 

rX2,     0<s<t, 
X(s)  =  <  X2  ,    t1<s<t2  (6) 

.  X3  ,    t2<s^t 

and  assuming  that  the  majority  of  the  crew  is  on  liberty 

over  the  weekend,  then  a  reasonable  intensity  function  for 

weekends  might  be 

'X2,     O^s^tj' 
X(s')  =  i  (7) 

,A3,   t1,<s<t 

where  X,  and  X3  are  identical  for  both  sets  of  equations. 

Thus  the  mean  value  function  is : 

t 
m(t)  =  [  X(s)ds  =  Xiti+X2 (t2-t1)+X3 (t  -t2)        (8) 


22 


and 


m(t')  =  j  X(s')ds  =  X2tx ,+A3(t-t1 ') 


(9) 


and  it  can  be  shown  that 

e-m(t)  [m(t)]n 

P{N(t)=n}  = ,   n>0 

n! 


(10) 


This  then  describes  a  non-homogeneous  Poisson  pro- 
cess.  The  values  for  each  X^  were  determined  by  comparing 
the  mean  value  function  m(t) ,  multiplied  by  4.9  gallons  of 
sewage/arrival,  to  the  expected  daily  and  weekly  sewage 
build-up.   The  details  of  this  are  contained  in  Appendix  A. 
The  final  values  used  for  each  Aj_  are  presented  in  Figure  3. 


Arrival/5  min 


X,=4.5 


1  A, =1.3 


A,  =  .5 


0500 


1600    2200    0500 


Time 


Figure  3.   Arrival  Rates  for  Simulation  One.  > 
The  simulation  involved  generating  Poisson  variates 
for  each  5  minute  increments  of  simulated  time,  using  the 
above  parameters  during  the  periods  indicated.   The  number 
of  arrivals  was  then  converted  to  a  quantity  of  sewage  X(t) 
according  to : 
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X(t)  =  N(t)  x  4.9  gal/arrival  (11) 

giving  the  necessary  random  input  for  X(t)  which  is  then 
applied  to  the  basic  model  as  described  in  section  1  above. 
3.   Input  Distribution  as  Described  by  Simulation  Two 

Although  the  reasoning  behind  Simulation  One  was 
considered  sound,  it  was  decided  that  a  second  simulation 
could  be  developed  through  a  different  approach.   Not  only 
would  it  take  full  advantage  of  all  known  data,  but  it  would 
provide  an  experiment  by  which  the  two  simulations  could  be 
weighed  against  one  another. 

In  order  to  get  the  best  approximation  to  the  true 
arrival  distribution  over  a  24  hour  period  on  board  ship,  a 
group  of  experienced  Naval  personnel  were  surveyed  for  their 
estimates  of  the  actual  arrival  pattern  (i.e.  hydraulic 
loading)  of  water  closets  and  urinals  over  a  24  hour  period. 
From  such  a  survey,  an  empirical  distribution  was  developed 
by  employing  a  technique  similar  to  that  used  to  get  job 
time  estimates  for  PERT  programming  (Program  Evaluation  and 
Review  Technique) .   The  details  of  this  method  of  estimation 
and  how  it  was  used  are  contained  in  Appendix  B.   The  re- 
sults are  shown  in  Figure  4a. 

Figure  4b  gives  a  similar  hydraulic  loading  pattern 
derived  for  a  commercial  cargo  ship  with  a  40  man  crew 
[Ref.  13,  p.  58].   Although  the  routine  on  a  commercial 
cargo  vessel  is  considerably  different  from  that  on  a 
military  vessel,  it  can  be  noted  that  "peak  loadings  are 
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double  the  average  loading  and  almost  five  times  the  minimum 
loading"  for  both  distributions,  and  that,  in  general,  they 
compare  favorably. 

Once  this  distribution  was  established,  it  was 
necessary  to  convert  it  to  realistic  quantities  of  sewage 
per  hour.   Again,  the  hydraulic  shipboard  data  obtained  in 
MEL  Report  346/6  4  [Ref.  19]  (which  is  the  basis  of  current 
Navy  design  parameters  of  30  g.m.p.d.)  was  used  to  estimate 
a  distribution  of  the  daily  per  capita  sewage  flow.   This 
data,  however,  was  revised  in  a  thesis  by  D.  B.  Campbell 
(a  fact  not  incorporated  in  design  specifications  listed  in 
Ref.  2,  p.  E-l) .   Combining  both  sets  of  information  it  was 
decided  that  daily  per  capita  sewage  flow  rates  could  best 
be  represented  by  a  beta  distribution  with  a  most  likely 
(mode)  value  of  14.8  gal/man/day,  a  maximum  of  20.0  g.m.p.d. 
and  a  minimum  of  11.3.'   As  shown  in  Appendix  C,  these  values 
were  used  to  estimate  the  parameters  of  a  beta  distribution 
[Ref.  17,  p.  196] . 

Observations  of  the  daily  sewage  generation  rate  per 
man  were  then  generated  from  this  beta  distribution  for  each 
day  of  the  simulated  week.   Once  a  daily  average  was  gen- 
erated it  was  converted  to  an  hourly  average  for  a  ship 
with  2  40  men.   (The  number  of  men  was  arbitrary  and  could 
be  read  in  by  the  programmer.)   Finally,  it  was  translated 
to  the  empirical  distribution  by  conversion  ratios  of  the 
hourly  relative  value  to  the  relative  mean.   For  example 
at  1200,  if  the  beta  realization  was  15  gal/man/day  then: 
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Relative  Value 

Hourly  Realization  = 

Relative  Mean 

Daily  Beta  Variate 
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x  No.  of  Men  (12) 


11.3    15 

Realization  =  x  —  x  240  =  232  gal 

7.3    24 

Such  values  were  the  hourly  realizations  of  X(t)  used 

as  the  input  to  the  basic  model  described  in  subsection  1. 
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IV.   RESULTS 

The  results  of  simulations  one  and  two  are  presented  on 
pages  32  through  35.   Figure  5a  and  6a  are  graphs  of  the 
average  maximum  level  of  sewage  held  as  a  function  of  pro- 
cessing rate  for  the  five  policies  tested  by  simulations  one 
and  two  respectively.   Figures  5b  and  6b  represent  the 
average  total  time  the  processor  must  run  in  a  week's  period 
(of  a  possible  168  hour  maximum) ,  again  as  a  function  of 
processing  rate  for  each  of  the  five  policies.   The  results 
agree  well  with  intuition. 

Looking  first  at  the  results  common  to  both  simulations, 
one  observes  that  the  maximum  average  sewage  level  has  an 
inverse  relationship  to  the  processing  rate.   The  percent 
savings  in  tank  size  decreases  with  higher  processing  rates, 
and,  for  rates  greater  than  210  gal/hr,  the  savings  appear 
to  be  negligible.   For  policies  in  which  processing  time 
can  vary,  the  total  weekly  processing  time  is  also  inversely 
proportional  to  the  processing  rate,  and,  over  the  range  of 
rates  considered,  the  function  is  nearly  linear. 

The  different  policies  considered  have  a  significant 
effect  on  sewage  level  and  processing  time.   Policy  two 
describes  a  limit  on  minimum  sewage  level  obtainable  for 
the  rates  considered.   The  curve  of  maximum  sewage  level 
increases  with  each  higher  policy  number  up  to  policy  five. 
Conversely,  for  average  processing  time  the  trend  is  just 
the  opposite  as  was  expected.   Policy  two  requires  the 
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maximum  processing  time  and  each  succeeding  policy  number 
requires  less.   For  the  policy  of  strictly  holding  sewage 
(policy  1) ,  the  average  sewage  accumulated  in  a  week's 
period  (assuming  the  ship  to  be  a  destroyer  with  a  comple- 
ment of  240  men)  was  21,100  gallons,  with  a  variance  of 
8.2  x  lO1*  gallons2  for  simulation  one,  and  21,080  gallons 
with  a  variance  of  23.4  x  10 5  gallons2  for  simulation  two. 
The  average  values  for  a  12  hour  holding  period  were  3,057 
gallons  for  simulation  one  and  2,700  gallons  for  simulation 
two. 

Certain  additional  observations  could  be  made  that  were 
peculiar  only  to  simulation  one.   First,  a  definite  change 
in  the  slope  of  the  curve  for  maximum  tank  level  vs  pro- 
cessing rate  occurs  at  150  gal/hr,  with  little  additional 
savings  in  tank  size  being  realized  beyond  that.   For  pro- 
cessing rates  above  150  gal/hr,  the  average  maximum  tank 
level  was  never  above  2,500  gallons.   The  variance/mean 
ratio  for  maximum  tank  level  is  on  the  order  of  from  2  to  10 , 
It  is  usually  highest  for  a  processing  rate  of  140  gal/hr, 
and  drops  off  immediately  thereafter.   For  policy  four  the 
ratio  is  much  lower  than  for  other  policies.   For  processing 
time  vs  processing  rate  the  variance/mean  ratio  is  on  the 
order  of  10~2  ,  and  policies  3  and  4  show  little  appreciable 
difference  in  time  savings . 

For  simulation  two,  it  can  be  seen  that,  with  policy 
two  as  a  reference,  the  curve  of  maximum  level  vs  processing 
rate  levels  off  more  rapidly  for  policies  3  and  4  beginning 
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at  approximately  160  gallon/hr.   The  curves  for  policies  3 
and  4  follow  the  same  contour  with  policy  3  having  a  maxi- 
mum sewage  level  on  the  order  of  1000  gallons  less  than 
policy  four  for  all  rates.   The  variance/mean  ratio  for 
maximum  sewage  level  is  on  the  order  of  102  to  10 3  at  140 
gal/hr  and  decreases  as  processing  rate  increases.   At 
higher  processing  rates  the  ratio  is  significantly  less  for 
policy  four  (approx.  20)  than  for  other  policies.   For  pro- 
cessing time  vs  processing  rate,  policies  three  and  four 
have  the  same  slope,  with  policy  four  showing  only  a  slight 
savings  in  time  of  about  5  hrs  for  the  week.   The  variance/ 
mean  ratio  for  either  policy  is  again  very  small,  on  the 
order  of  0.5. 

A  comparison  of  the  two  simulations  reveals  that,  al- 
though the  average  results  of  simulation  one  are  on  the 

i 

same  order  as  those  of  simulation  two,  the  variance  of  the 

results  is  much  smaller  for  the  former.   This  is  not  supris- 
ing,  since  simulation  one  is  based  on  the  Poisson  process 
in  which  the  variance  of  the  asymptoptic  distribution  of 
N(t)  will  equal  the  mean.   In  the  case  of  a  non-homogeneous 
process  in  which  the  parameters  do  not  differ  greatly,  it 
can  also  be  expected  that  as  the  number  of  replicates  in- 
creases the  sample  variance  will  approach  some  value  on  the 
same  order  of  magnitude  as  the  mean.   For  simulation  two, 
the  realizations  of  the  daily  average  sewage  rate  per  man 
can  vary  from  11.3  gal/man/day  to  20.0  gal/man/day  which  is 
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then  imposed  on  an  empirical  distribution  with  even  greater 
variance.   The  result  is  a  process  with  reasonably  large 
variance;  certainly  much  greater  than  that  resulting  from 
simulation  one. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.   INTERPRETATION  OF  RESULTS 

In  order  to  prescribe  the  most  feasible  combination  of 
processing  rates,  policies  and  holding  tank  size,  more 
information  regarding  the  cost  function  of  each  particular 
parameter  would  have  to  be  developed.   If  costs  were  as- 
sumed to  follow  some  increasing  linear  function  of  pro- 
cessing rate,  holding  tank  size  and  processing  time,  then 
the  policy  and  rate  which  would  minimize  a  measure  of  ef- 
fectiveness taking  all  such  costs  into  account  would  be 
expected  to  be  most  feasible. 

Based  on  the  results  of  Chapter  IV  and  with  the  under- 
standing that  cost  functions  must  be  applied,  it  appears 
that  a  system  consisting  of  a  holding  tank  of  2,700  gallon 
capacity  for  sewage  combined  with  a  processing  unit  capable 
of  handling  160  to  170  gallons  of  sewage  per  hour  would  be 
most  beneficial.  Further  such  a  system  should  be  operated 
under  a  processing  policy  similar  to  that  of  policy  3. 

It  was  contended  in  Chapter  I  that  the  Navy's  design 
parameter  of  30  gallons/man/day  was  over  estimated.   The 
simulations,  which  were  developed  using  the  best  data  that 
could  be  obtained,  bear  out  this  conclusion.   More  specifi- 
cally, if  this  current  value  were  used  to  construct  a  12 
hour  holding  tank  for  a  destroyer  of  2  40  men,  the  minimum 
tank  size  would  be  3600  gallons.   The  values  obtained 
through  simulation,  which  were  measured  during  the  period 
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when  traffic  flow  was  the  greatest  (0500  to  1700),  were, 
again,  3,050  gallons  for  simulation  one,  and  2,700  gallons 
for  simulation  two.   Even  with  the  addition  of  a  200  gallon 
safety  factor,  it  appears  evident  that  a  savings  of  between 
400  to  700  gallons  could  be  realized  for  each  12  hour 
holding  tank  installed.   Even  greater  savings  would  be 
expected  for  larger  tanks. 

The  capacity  of  a  holding  tank  based  on  the  quantity  of 
diluent  (flushing  liquid)  used  can  also  be  computed  using 
an  equation  presented  in  a  report  by  the  Maritime  Research 
Information  Service  [Ref.  13].   Namely: 

C  =  1.25  x  M  x  (B+D)  (13) 

where 

C  =  tank  capacity 

M  =  man  days  of  storage 

B  =  daily  per  capita  body  waste  =  0.437  gallons 

D  =  diluent  in  gallons/man/day. 

If  the  value  of  diluent  is  taken  to  be  15  gal/man/day, 

the  unmodified  mean  of  the  data  obtained  by  MEL  Report 

346/64,  then  the  recommended  holding  tank  size  would  be 

2,320  gallons.   The  values  obtained  by  simulation  would 

thus  appear  conservative;  however,  it  is  important  to 

realize  that  they  also  account  for  peak  arrival  periods. 

B.   RECOMMENDATIONS 

The  simulations  described  in  this  thesis  offer  sufficient 
merit  for  use  in  the  design  of  sewage  treatment  systems.   It 
is  recommended  that  a  study  be  made  to  obtain  further 
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statistics  on  the  arrival  pattern  of  head  facilities  aboard 
a  Naval  ship.   Such  new  data  could  then  replace  the  empirical 
distribution  developed  for  simulation  two,  or  be  used  to 
revise  the  intensity  function  of  simulation  one.   This 
should  lead  to  even  more  meaningful  results  than  those 
presented . 

It  is  further  recommended  that  design  parameter  for 
daily  per  capita  sewage  generation  be  reviewed.   If  it  is 
decided  that  it  is  in  fact  overestimated,  then  it  is  recom- 
mended that  12  hour  holding  tanks  of  the  size  presented  in 
this  report  be  tested  on  board  destroyer-sized  ships,  and 
that  the  simulations  be  employed  to  develop  holding  tanks 
for  larger  ships  as  well. 
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APPENDIX  A 
CALCULATION  OF  ARRIVAL  PARAMETERS 

As  pointed  out  in  the  body  of  the  thesis,  the  mean  value 
function  was  used  to  estimate  the  best  values  of  X^  (i= 
1,2,3) .   This  was  done  by  computing  average  daily  and  weekly 
sewage  generation  first  by  using  the  average  values  pre- 
sented in  MEL  Report  346/64  [Ref.  19],  and  second  by  use  of 
policy  one  and  simulation  two.   The  values  obtained  from 
weekly  observations  of  sewage  build-up  as  estimated  by 
simulation  two  were  used  since  1)  they  were  considered  more 
representive  and  2)  they  would  also  serve  as  a  means  of 
normalizing  both  simulations  so  that  over  long  periods  of 
time  the  two  simulations  would  give  similar  results  for 
strict  holding.   Average  daily  and  weekly  sewage  generation 
was  then  divided  by  4.9  gal/arrival  to  determine  the  number 
of  arrivals/period. 

The  results  were : 

Weekly  average  =  21,000  gal/week/4.9  gal/arrival 

=  4260  arrivals/week 
similarly : 

daily  average  =  730  arrivals/day. 
These  were  then  equated  to  the  mean  value  function  which 
was 

M(t)  =  5(m(t))  +  2  (m(t'))  =  5 [X 1 (t , ) +  X2 (t2 -tx ) +  X 3 (t-t2 ) ] 
+  2[X2(t1  ')+X3  (t-t/)]  (A-l) 
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for  the  weekly  value,  and 

m(t)  =  \1(t1)+\2(t2-t1)+\3(t-t2)  (a-2) 

for  the  daily  value. 

Substituting  in  the  values  for  tj_  where  t  =  24  hours  = 
28  8  five  minute  increments  and  setting  them  equal  to  the 
weekly  and  daily  sewage  build-ups,  one  gets: 

weekly:   660A!  +  768A2  +  588A3  =  4260  (A-3) 

daily:    132AX  +  72A2  +  84A3  =  730  (A-4) 

where:    t1    =  132,   t2  =  204,   t  =  288,  t *  =  132. 

In  order  to  solve  the  above  simultaneous  equations  a 
third  equation  was  needed.   This  was  obtained  from  weekends 
and  computing  a  third  average  build-up  for  that  period  to 
obtain : 

2  m(t')  =  2,100/4.9  =  420  arrivals/weekend 
or 

264A2  +  312A3  =  420.  (A-5) 

Solving  (A-3)  through  (A-5)  simultaneously  one  obtains 
the  final  estimates  of  A-  (i  =  1,2,3)  namely: 

Aj  =  4.5  arrivals/5  min 

A2  =  1.0  arrivals/5  min 

A 3  =   .5  arrivals/5  min 

which  agree  well  with  intuition. 


40 


APPENDIX  B 
METHOD  OF  ESTIMATION  FOR  AN  EMPIRICAL  DISTRIBUTION 

In  order  to  obtain  the  hydraulic  loading  pattern  of 
water  closets  and  urinals  aboard  Naval  ships,  it  was  nec- 
essary to  resort  to  a  subjective  estimation  technique. 
The  technique  used  is  the  same  as  that  for  obtaining  job 
time  estimates  for  PERT  scheduling  computations  [Ref.  14]. 
Such  a  technique  is  employed  when  one  is  unable  to  do  any 
statistical  sampling. 

It  was,  therefore,  decided  to  survey  a  cross-section 
of  experienced  Naval  personnel,  both  officers  and  enlisted, 
for  their  subjective  estimates  of  the  relative  usage  of 
head  facilities  for  each  hour  of  a  24  hour  period.   Their 
estimates  were  made  on  a  relative  scale  of  0  to  20 ,  (0  in- 
dicating the  facility  was  not  used  at  all,  20  indicating 
it  was  used  to  capacity) .   Each  person  was  to  give  a  high, 
most  likely,  and  low  estimate,  independent  of  estimates 
made  by  others . 

It  has  been  "historically  accepted"  that  job  duration 
times  described  by  such  estimates  are  beta  distributed  over 
a  finite  interval  given  by  the  high  and  low  estimate.   The 
same  distribtuion  can  be  assumed  to  apply  to  estimates  of 
the  arrival  pattern.   The  expected  value  of  such  a  distri- 
bution can  then  be  approximated  by: 


Ue  =  (a+4m+b)/6  (B-l) 
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and  the  variance  by; 

Va  =  I(b-a)/6]2  (B-2) 

where : 

a  =  low  (optimistic)  estimate 

m  =  most  likely  (mode)  estimate 

b  =  high  (pessimistic)  estimate. 
Using  the  expected  value,  Ue ,    as  the  best  estimate  of  each 
individual .   All  such  values  were  then  averaged  to  give  the 
overall  best  estimate  for  each  hour  of  the  24  hour  period. 
The  results  are  the  empirical  distribution  plotted  in  the 
body  of  the  thesis  (Figure  4a) . 
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APPENDIX  C 
ESTIMATION  OF  BETA  PARAMETERS 

A  beta  distribution  that  is  transformed  to  some  interval 
other  than  (0,1)  is  given  by: 


r(A+B) 

(b-a)1^"6  (T-a)A~1  (b-T)5"1,  b>T>a 

r(A)T(B) 

(C-l) 
0  elsewhere 


f(T)  = 

where  T  is  the  transformed  beta  variable  given  by: 

T  =  a+(b-a)X  (C-2) 

X  being  a  (0,1)  beta  variable  and  a,b  being  the  lower  and 
upper  bound  of  the  transformed  variable.   If  the  high  (b) , 
low  (a)  and  most  likely  (m)  values  of  the  distribution  can 
be  determined,  then  the  transformed  beta  parameters  can  be 
related  to  (0,1)  beta  parameters  according  to: 

A  =  46E(X)  [E(X)-E(X)2  -  ^}  (C-3) 

A-AE(X) 

B  =  (C-4) 

E(X) 

where 


4m+b-5a 

E(X)  =  (C-5) 

6(b-a) 

Using  the  values  obtained  from  MEL  Report  346/64  i.e.: 


a  =  11.3  gal/man/day 
m  =  14.8  gal/man/day 
b  =  20.0  gal/man  day. 
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Then  the  values  for  the  parameters  can  be  solved  for,  using 
the  above  equations  to  obtain: 

A  =  4.3       B  *  5.7 
These  values  are  then  used  to  generate  a  beta  variable,  X, 
by  first  generating  two  independent  gamma  variates  Y^GCl^A) 
and  Y2(1,B)  and  computing  the  ratio 


X  =  (C-6) 

Y1+Y2 

It  then  remains  only  to  transform  X  to  a  realization  of  T 

using  equation  (C-2) . 


44 


APPENDIX  D 
USE  OF  SIMULATION  ONE  AND  TWO 

The  two  simulations  were  programmed  in  Fortran  IV  for 
use  on  the  IBM  360/6  7  computer.   Both  simulations  allow  the 
option  of  plotting  the  graph  of  the  sewage  build-up  (in- 
ventory) for  each  hour  of  the  week.   (Examples  of  such 
plots  are  given  under  computer  output.)   If  one  uses  the 
option,  however,  he  must  insure  that  the  appropriate  JCL 
cards  are  used  for  his  particular  hardware. 

Only  one  data  card  is  needed  for  either  simulation. 
All  that  must  be  provided  are  the  parameters  to  be  tested, 
an  indication  of  print  routine  to  be  used  and  a  flag  for 
calling  the  plot  routine.   For  each  simulation  the  format 
of  the  data  card  is  as  follows : 

a)   Simulation  One.   The  following  variables  must 

be  supplied: 

Variable     Cols  on  Card        Description 

As  described  in  the  body 
of  the  thesis,  these  are 
the  arrival  rates/5  min 
that  are  used  for  the 
Poisson  generator. 

The  initial  processing 
rate  seed,  in  gal/5  min 
processed.   It  is  in- 
creased by  10  gal/hr  as 
soon  as  it  is  read  in. 

NPLAN  19-20       A  two  digit  integer 

specifying  the  processing 
plan  to  be  tested. 


Plan   1 

0-5 

Plan   2 

6-10 

Plan    3 

11-15 

Alpha   1 

16-18 
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NOPRIN  21-22       A  two  digit  integer  to 

indicate  which  of  four 
print  routines  is  to  be 
used.   (See  below  for 
discription  of  these 
routines . ) 

NO  DRAW         23-24       A  flag  to  indicate  if 

plot  routine  is  desired 
(any  2  digit  number)  or 
not  (use  zero)  . 

ZSTOP,ZSTART     25-20        These  are  the  level 

30-3  4       parameters  for  policy  4 
in  F5 . 1  format . 

TSTOP  34-36       A  3  digit  integer  repre- 

senting the  stopping 
time  for  policy  5. 

For  simulation  two  the  same  values  are  read  in  with  the 
exception  of  the  Poisson  parameters,  instead  the  number  of 
men  representing  the  size  of  ship  is  read  in  as  a  3  digit 
integer. 

For  the  print  routines  it  should  be  noted  that  if  full 
printout  is  used  that  this  can  add  considerably  to  the  run- 
ning time  of  the  program.   Depending  on  the  read  in  for 
NOPRIN ,  the  routines  are : 

NOPRIN  ROUTINE 

1  Print  only  the  final  statis- 
tics for  ZMAX  matrix  and  the 
processing  time  matrix. 

2  Print  above  plus  the  times 

that  the  processor  stop  and 
restarted  during  the  week. 
Print  also  the  processing 
rate ,  run  number  and  maximum 
tank  level  on  each  run. 

3  Print  the  above  plus  the 

entire  Z(t)  matrix  (also  the 
X(t)  matrix  for  simulation  2) 
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4  Print  the  above  plus  the 

X(t)  matrix  and  arrivals 
matrix.   (Simulation  One) 

In  conclusion,  it  should  be  pointed  out  that  simulation 
one  requires  approximately  2  000  random  numbers  to  be  gen- 
erated per  replication;  it  requires  considerably  more  com- 
puter time  than  simulation  two.   This  is  inpart  compensated 
for,  however,  by  its  smaller  variance,  so  that  fewer  runs 
are  needed. 
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COMPUTER  OUTPUT 


(Sample    only) 
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COMPUTER    PROGRAM 


*  TANK    SIMULATION    ONE  * 

*JLf    %f--    «JL-    <JU   •Xr    -Xf    \Xr    -Xr    «L>    -»*^    -J«r    «Lr     «X-    «.«.,    »<,     ^t,  <JLi    OL-    •*'  •    -Xj    -X,    tXr    -X-   OU     -X-    >!*     ai>    -Xf     *A«    "A"    J-     *>.-    <JL>    J-    -A-   >A»    Jrf    Jf    ■•*-    -J*    «!«■    »Xr    sAr    "«V    *JLr   J*   aJL>    Jt*    Xr    -X-r    -Xr     -%\r    -X-    *J#    -J--    <JL>     vV    ■"■'•'    «J»    — L" 

Jfi   -v    -,-    -v   ^i    A^.   <Y>   T    -v-    -,-    —    ^  ~    ^    ,,.   ^.    ^.   0>   ,,,    ^„    ,,^    ,, ,,    ,,,   -p     -,-    •,-    -,-     -,-    -,-   *p    -,-    -v    -,-.    n-    -v*  -p.    ~,-    -^    -t>    -,-   ',-    *,-    -,~    -p    -,-    ^    -,.,,,,.    .^    „,.,,-.   .,n    ^p    »^    ^    ,*.    J^    n„    ^ 

THIS    IS    THE    FIRST    OF    TWO    SIMULATIONS    DEVELOPED    BY    J.O.MINER, 
JR.     TO    DESCRIBE    THE    PROCESS    OF    GENERATING    AND    TREATING 
SEWAGE    ON    BOARD    NAVAL    SHIPS.     IT    EMPLOYS    A    NON-HOMOGENEOUS 
POISSON    PROCESS    TO    DESRIBE    THE    ARRIVAL    OF    SEWAGE    TO    A    HOLD- 
ING   TANK    OVER    THE    COURSE    OF    THE    DAY.     IT    THEN    TESTS    FIVE    PRO- 
CESSING   POLICIES    FOR    ANY    PROCESSING    RATE    DESIRED.    THE    SIZE 
OF    THE    SHIP    IS    ARBITRARY    AND    IS    READ    IN    BY    THE    PROGRAMMER. 

j,  j.  O,  x  j,  g,  j,  j,  o,  o.  v1.  j.  j,  f  ^  *i,  j,  o,  v-  a,  j,  ^  j,  j,  o-  J,  v,  J-  x  J-  J.  J,  .1-  f  J/  »u  v.  J-  >>,  o,  j.  ^.  j.  j,  j,  vi,  j,  x  .1,  j,  j.  j,  j,  j.  j,  x  ,1,  a.  x  x 

.  ?,■■   *y*.   *,*.  ?{±  ^-.  *,-.  ^^  »^«.  *f%  •»(*.  *vs  ^r»  *C'  *f  *S*  *t*  *?■  "*i*  ^*  *r"  't*  *V*  *t*   *t*  T**  "T*  *>   *P  *TT  *P"  ^  *!*"  *P  *T*  •!*  T^  *•*"  *^  *l*  T*  T*  *1*  *P  *F"  *"P  TT  *1%  *T*  ""<*  "T*  ^P  HP  *¥*  *"P  T"  HP"  *1P  *Y*  *¥*•  T* 

*  * 

THE    FOLLOWING    STATEMENTS    INITIALIZE    THE    PROGRAM    AND    READ    IN 
THE    VARIOUS    PARAMETERS    TO    BE    TESTED, INCLUDING    THE    PARAMETERS 
FOR    THE    POISSON    DISTRIBUTION 

COMMON    I,  J,ITIME,NOPRIN 

DIMENSION    Z(7,288) ,X(7,288) ,NUM(7,288) 

DIMENSION    XORD( 168 ) , YORD ( 168 ) , ZMAX( 10,20) , TPROSS ( 10 , 20 
REAL*8    TITLE(12)/"TANKSIM1",11*"  '/ 

REAL    LABEL/'  ■/ 

INTEGER    TSTGP 

RE  AD  (5, 100)    PLAM1,PLAM2,PLAM3,ALPHA1,NPLANJ,N0PRIN, 
CN0DRAW,ZST0P,Z START, TSTOP 
100    F0RMAT(3F5.3,F3.1,3I2,2F5.1,I3) 
IX=4321571 

a,  a* 

RUNS    ARE    MADE    FOR    NINE    PROCESSING    RATE    IN    INCREMENTS    OF    10 

GALLON/HOUR 

t~  £ 

DO    6000    M=l,9 
ALPHAl=ALPHAl+.833 
DO    5000    L=l,10 
IFINOPRIN.LT. 2)    GO    TO    155 
WRITE(6,150)L 
150    FORMAT  PI"  ,20X,"RUN    NO.       •  ,12,//) 

VARIOUS    ARRAYS    ARE    ZEROED    FOR    EACH    RUN 

155    ZMAX(M,L)=0.0 

ALPHA=0.0 

ITIME=0 

DO    30    1=1,7 

DO    20    J=l,288 

Z(I,J)=0.0 

X{ I ,J)     =    0.0 

NUM( I , J)    =    0 
20    CONTINUE 
30    CONTINUE 

IX    =       IX    +    342 

A    POISSON    ARRIVAL    IS    THEN    GENERATED    FOR    EACH    FIVE    MINUTE 
INCREMENT    OF    THE    SIMULATED    WEEK.     DIFFERENT    PARAMETERS    ARE 
USED    ACCORDING    TO    THE    TIME    OF    DAY    OR    WEEK. 

DO    600    1=1,5 

PLAM    =    PLAM1 

DO    500    K=l,132 

CALL  P01SSN(PLAM,IX,N) 
500    NUM(I,K)=N 

PLAM    =PLAM2 

DO    510    K=133,204 

CALL  P01SSN(PLAM, IX,N) 
510    NUM( I ,K)=N 


52 


PLAM  =  PLAM3 

DO  520  K=205,288 

CALL  P01SSM(PLAM,IX,N) 
520  NUM( I,K)=N 
600  CONTINUE 

DO  750  1=6,7 

PLAM=  PLAM2 

DO  700  K=l,204 

CALL  P01SSN(PLAM,IX,N) 
700  NUM(I,K)=N 

PLAM=PLAM3 

DO  720  K=205,288 

CALL  P01SSN(PLAM, IX, N) 
720  NUM(I,K)=N 
750  CONTINUE 

DO  1500  1=1,7 

THF  ARRIVALS  ARE  THEN  CONVERTED  TO  A  QUANTITY  OF  SEWAGE  AND 
THE  SEWAGE  PROCESSED  IN  ACCORDANCE  WITH  ONE  OF  THE  FIVE  POL- 
ICIES. 

IF( (NPLAN.EQ.l) .OR. (NPLAN.EQ.4) )  GO  TO  90 

ALPHA=ALPHA1 
90  DO  1000  J=l,288 

X(I  ,J)=NUM( I, J)*4.94 

IF((  J.EQ.D.AND.  (I.EQ.l)  )G0  TO  855 

IF((  J.EQ.D.AND.  (I.GE.2)  )G0  TO  860 

Z(I,J)=  Z( I,J-1)+X( I,J J-ALPHA 

GO  TO  880 
855  IF  (NPLAN.EQ.4)  ALPHA=ALPHA1 

Z( I, J)=X( I,J)-ALPHA 

GO  TO  880 
860  Z(I,J)=  Z(I-1  ,288)+X( I,J)-ALPHA 
880  IF(Z( I,J).LT.O.O)Z( I, J  )  =  0.0 

IF(ZMAX(M,L) .LT.Zt I, J) )  ZMAX ( M, L )=Z( I , J ) 

TANK=Z(I, J) 

IF(NPLAN.EQ.l)  GO  TO  1000 

IF(NPLAN.EQ.2)  GO  TO  1000 

IF(NPLAN.EQ.3)CALL  PL AN3 ( TANK , AL PHA1 , ALPHA , £900) 

IF (NPLAN.EQ.4) CALL  PLAN4(TANK, AL PHA1 , ALPHA , ZSTOP , 
CZSTART,£900) 

IF(NPLAN.EQ.5)  CALL  PLAN5( ALPHA1 , ALPHA , TSTOP , £900 ) 
900  CONTINUE 
1000  CONTINUE 
1500  CONTINUE 

THE  WEEKLY  PROCESSING  TIME  IS  THEN  COMPUTED 

IF(NPLAN.EQ.l)  ITIME=2016 
IF(NPLAN.EQ.2)  ITIME=0 
IF(NPLAN.EQ.5)  IT  I ME= ( 288-TSTOP ) *7 
TPROSS(M,L) =168.-1  TIME/1 2. 

IF  DESIRED,  THE  WEEKLY  SEWAGE  INVENTORY  MAY  BE  PLOTTED. 

IF(NODRAW.EQ.O)  GO  TO  2900 

INDEX=0 

DO  2000  1=1,7 

DO  1900  J=12,288,12 

INDEX=INDEX+1 

YORD( INDEX)=Z( I, J) 
1900  CONTINUE 
2000  CONTINUE 

DO  2100  J=l,168 
2100  XORD( J  )=J 

CALL  DR AW ( 168, XORD,YORD, 0,0, LAB  EL, TITLE, 20. 0,300. 0,0,0 
1LAST) 

REALIZATION  OF  EITHER  SEWAGE  OR  ARRIVALS  FOR  EACH  FIVE  MIN-" 
UTE  INTERVAL  MAY  THEN  BE  PRINTED  AS  DESIRED.  Z  VALUES  ARE 
CUMMULATIVE  HOURLY 
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2900  IF(N0PRIN.LT.2)  GO" TO  5000 

WRITE (6, 3500 )ZMAX(M,L)  , ALPHA  1 , TPRO SS ( M , L ) 
3500  FORMATCO'  ,20X,«  ZMAX=«  ,F8.1,10X,  'PROCESSING  RATE  =  ', 

CF6.1,10X, 'PROCESSING  TIME  =  «,F7.2,//) 
2950  IF  (N0PRIN.LT.3)  GO  TO  5000 

WRITE (6,3000) ( ( Z( I , J ) , 1=1, 7 ) , J=12, 288, 12) 
3000  FORMATCO   Z  MATRI X • // ( IX, 7 ( F8. 1 ,4X) ) ) 
4000  IF(N0PRIN.LT.4)  GO  TO  5000 

WRITE (6, 3100) UXU, J), 1  =  1, 7), J  =  l, 288) 
3100  F0RMATC1   X  MATRI X • // ( IX, 7 ( F5. 2 ,3X) ) ) 

WRITE (6,3200) ( ( NUM ( I, J), 1=1, 7), J=l, 2 88) 
3200  FORMAT  CI   ARRIVALS  MATRIX • //( IX ,7 ( 12, 3X )  )  ) 
5000  CONTINUE 
6000  CONTINUE 

STATISTICS  ARE  THEN  COMPUTED  FOR  MAXIMUM  SEWAGE  LEVEL  AND 
FOR  PROCESSING  TIME.  THE  FINAL  RESULTS  ARE  THEN  PRINTED 

DIMENSION  ZSUM(20) ,ZAVE(20) ,ZVSUM{ 20) ,TSUM(20) ,TAVE(20 
C),TVSUM(20),ZVAR(20) ,TVAR(20) 

DO  7200  J=1,M 

ZSUM( J )  =  0. 

ZVSUM( J)=0. 

TSUM( J)=0. 

TVSUM( J)=0. 

DO  7000  1=1, L 

ZSUM( J)=ZSUM( J)+ZMAX(J,I  ) 

TSUMU  )=TSUM(J)+TPROSS(J,I  ) 
7000  CONTINUE 

ZAVE(J)=ZSUM(  J)/L 

TAVE( J)=TSUM(J)/L 

DO  7100  1=1, L 

ZVSUM(  J)=ZVSUM( J)+( (ZMAX( J.I)-ZAVE(J) )**2) 

TVSUM( J)=TVSUM( J )+ ( ( TPROSS ( J , I )-TAVE( J) )**2) 

ZVAR( J)=ZVSUM( J)/L 

TVAR(J)=TVSUM'(J)/L 
7100  CONTINUE 
7200  CONTINUE 

WRITE(6,8000)  NPLAN 
8000  FORMAT  CI     ZMAX  MATRIX  FOR  POLICY  NO.  ',12//) 

DO  8200  J=1,L 

WRITE(6,8100)( ZMAX( I TJ ) ,I=1,M) 
8100  FORMAT ( IX, 10 ( F8 .2, 3X  )  ) 
8200  CONTINUE 

WRITE (6, 8500) (ZAVE( J) , J=1,M) 
8500  FORMATCO     AVERAGE  ZMAX  •  /  /  (  3X  ,  10  (  F8.  2  ,  2X  )  )  ) 

WRITE (6, 8600) ( ZVAR( J) , J=1,M) 
8600  FORMATCO     ZMAX  VARIANCE  « // (  3X,  10(  F  10.  2  ,  IX  )  )  ) 

WRITE(6,9000)  NPLAN 
9000  FORMAT  CI    PROCESSING  TIME  MATRIX  FOR  POLICY  NO.  «, 
CI2//  ) 

DO  9200  J=1,L 

WRITE ( 6, 8150) (TPROSS (I, J), I =1,M) 
8150  FORMAT(1X,10(F7.2,3X)  ) 
9200  CONTINUE 

WRITE (6, 9300) ( TA VE ( J ) , J=l , M ) 
9300  FORMAT(»0     AVERAGE  T IME ' / / ( 3X , 10 ( F6. 2 , 2X ) ) ) 

WRITE (6,9400) (TVAR( J) , J=1,M) 
9400  FORMATCO     TIME  VARIANCE  •  //  (  3X,  10(  F  10.  2  ,  IX)  )  ) 

STOP 

END 

THE  FOLLOWING  SUBROUTINE  GENERATES  POISSON  VARIATES  FOR  EACH 
FIVE  MINUTE  INCREMENT.  A  TABLE  LOCK-UP  METHOD  WAS  USED 
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SUBROUTINE    POISSN ( PL AM , I X,N ) 

COMMON     I,J,ITIME,NOPRIN 

INTEGER    N 

A=EXP(-PLAM) 

B=A*PLAM+A 

C=( A*PLAM**2)/2+B 

D=(A*PLAM**3)/6+C 

E=( A*PLAM**4)/24+D 

F=(A*PLAM**5)/120+E 

G=(A*PLAM**6)/720+F 

H=( A*PLAM**7)/5040.+G 

P=< A*PLAM**8) /4032  0.+H 

5  CALL    RANDUUX,  IY,YFL) 
IX=IY 

R=    YFL 

IF(R.LE.A)  GO  TO  6 
IE  (R.LE.B)  GO  TO  7 
IF  (R.LE.C)  GO  TO  8 
IF  (R.LE.D)  GO  TO  9 
IF  (R.LE.E)  GO  TO  10 
IF  (R.LE.F)  GO  TO  11 
IF  (R.LE.G)  GO  TO  12 
IF  (R.LE.H)  GO  TO  13 
IF  (R.LE.P)  GO  TO  15 
N=9 
GO  TO  14 

6  N=0 

GO  TO  14 

7  N=l 

GO  TO  14 

8  N=2 

GO  TO  14 

9  N=3 

GO  TO  14 

10  N=4 

GO  TO  14 

11  N=5 

GO  TO  14 

12  N=6 

GO  TO  14 

13  N=7 

GO  TO  14 
15  N=8 

14  RETURN 
END 

SUBROUTINE  PL AN3 (TANK , ALFA1 , RATE ,* ) 

THIS  POLICY  PROCESSES  SEWAGE  UNTIL  1600  AFTER  WHICH  TIME  IT 
THEN  EXAMINES  THE  LEVEL  OF  THE  TANK  AND  IF  BELOW  50  GAL. IT 
STOPS  PROCESSING  UNTIL  0500  THE  FOLLOWING  DAY.  ON  WEEKENDS 
PROCESSING  STOPS  ANYTIME  THE  LEVEL  IN  THE  TANK  REACHES  50  GA 

COMMON  I,  J,ITIME,NOPRIN  , 

IF((  I  .  EQ.l)  .AND.U.EQ.  1)  )  IFLAG=0  ' 

IF(N0PRIN.LT.2)  GO  TO  6 
IF(( I .EQ.l). AND. ( J.EQ. 1) )  WRITE (6,5) 

5  FORMATCO     PROCESSING  POLICY   NO. 3  WAS  USED  ON  THIS 
CRUN« ,//) 

6  IF( ( J. EQ.l). AND. ( IFLAG.EQ. 1) )  I T IME=ITI ME+288-JST0P 
IF(J.EQ.l)  IFLAG=0 

IFU.LE.5)  GO  TO  10 
IF(TANK.LE.50. )  GO  TO  30 

IF( (RATE. LT. 0.01 ) .AND. (TANK. GE  .  500.) )  GO  TO  20 
GO  TO  40 
10  IFU.LE.132)  GO  TO  40 

IF(TANK.LE.50. )  GO  TO  30 
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GO  TO  40 
20  RATE=ALFA1 

IF(IFLAG.EQ.O)  GO  TO  40 
ITIME=ITIME+J-JSTOP 
IF(N0PRIN.LT.2)  GO  TO  26 
WRITE(6,25)I  ,J 

25  FORMATC     CONTINUED  PROCESSING  AT  T=«  ,  II  ,  •  -•  ,  13  ,  //  ) 

26  IFLAG=0 
GO  TO  40 

30  RATE=0.0 

IFdFLAG.EQ.DGO    TO    40 
JSTOP=J 

IF(N0PRIN.LT.2)  GO  TO  36 
WRITE(6,35)I ,J 

35  FORMATCO      STOPPED  PROCESSING  AT  T=  ■  *  II  *  •-  ■  ,  13,  //  ) 

36  IFLAG=1 

40  IF( ( I.EQ.7).AND. (J.EQ.288) )  GO  TO  45 

GO  TO  SO 
45  IF( IFLAG.EQ.l)  I TI ME=I TIME+268- JSTOP 
50  RETURN1 

END 

SUBROUTINE  PL AN4 (TANK , ALFA1 , RATE , STOP, START,*) 

THIS  IS  A  LEVEL  CONTROL  POLICY  WHICH  CONTINUOUSLY  EXAMINES 
TG&  L&OEL  OF  SEWAGE  IN  THE  TANK.  IF  BELOW  SOME  SET  LEVEL 
(ZSTOP)  THEN  PROCESSING  STOPS  UNTIL  THE  SEWAGE  LEVEL  REACHES 
SOME  UPPER  VALUE  (ZSTART)  AT  WHICH  TIME  PROCESSING  CONTINUES 

COMMON  I ,J,ITIME,NOPRIN 

IF(  (  I.EQ.D.AND.  (  J.EQ.  1)  )  IFLAG  =  0 

IF(N0PRIN.LT.2)  GO  TO  6 

IF((  I  .EQ.D.AND.  (J.EQ.D)  WRITE(6,5) 

5  FORMATCO     PROCESSING  POLICY   NO. 4  WAS  USED  ON  THIS 
CRUN' ,//) 

6  IF(  (J.  EQ.D.AND.  (  IFLAG.EQ.l)  )  I TIME=I TI ME+288-JSTOP 
IF(J.EQ.l)  JSTOP=0 

IF(TANK.GT.STOP)  GO  TO  20 

RATE=0.0 

IFdFLAG.EQ.DGO.   TO    40 

JSTOP=J 

IF(N0PRIN.LT.2)  GO  TO  11 

WRITE(6,10)I,J 

10  FORMATCO      STOPPED  PROCESSING  AT  T=  •  ,  II  ,  ■-■  ,  13,  //  ) 

11  IFLAG=1 
GO  TO  40 

20  IF(TANK.LE. START)  GO  TO  40 
RATE=ALFA1 

IF( IFLAG.EQ.O)  GO  TO  40 
ITIME=ITIME+J-JSTOP 
IF(N0PRIN.LT.2)  GO  TO  31 
WRITE(6,30) I ,J 

30  FORMATC     CONTINUED  PROCESSING  AT  T=«  ,  1 1  ,  ■-•  ,  I  3  ,  //  ) 

31  IFLAG=0 

40  IF( (I . EQ. 7) .AND. (J.EQ.288) )  GO  TO  45 

GO  TO  50 
45  IF(IFLAG.EQ.l)  I T IME= ITIME+288- JSTOP 
50  RETURN1 

END 

SUBROUTINE  PL AN5 ( ALFA1 , RATE ,TIME,* ) 

THIS  IS  A  TIME  CONTROL  POLICY.  IT  ALLOWS  PROCESSING  FOR  A 
FIXED  PORTION  OF  THE  DAY  ONLY 

COMMON  I,J,ITIME,NOPRIN 
INTEGER  TIME 
IF(N0PRIN.LT.2)  GO  TO  6 
IF((  I  .EQ.D.AND.  (J.EQ.D  )  WRITE(6,5) 
5  FORMATCO     PROCESSING  POLICY   NO. 5  WAS  JSED  ON  THIS 
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CRUN* ,//) 
6  IF(J.LE.TIME)  GO  TO  10 
RATE=0.0 
10  RETURN1 
END 


SUBROUTINE  RANDU ( I  X, IY , YFL ) 
# 

THIS  ROUTINE  GENERATES  RANDOM  NUMBERS  FOR  USE  IN  THE  MAIN 
PROGRAM 

COMMON  I ,J,ITIME,NOPRIN 
IY=--!X*65  539 
IF(IY)5,6,6 

5  IY=I Y+2147483647+1 

6  YFL=IY 

Y FL=YFL*. 465661 3 E-9 

RETURN 

END 


Note:   Comment  cards  denoted  by  and  including  *  should  be 
removed  prior  to  running  programs . 
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*  TANK    SIMULATION    TWO  * 

THIS    SIMULATION    WAS    DEVELOPED    BY    JOHN    0    MINER, JR    TO    DESCRIBE 
THE    PROCESS    OF    GENERATING    AND    PROCESSING    SEWAGE    ABOARD    A 
NAVAL    SHIP.     IT    ASSUMES    THAT    THE    DAILY    AVERAGE    SEWAGE    GENERA- 
TION   RATE    PER    CAPITA    IS    BETA    DISTRIBUTED    OVER    THE     INTERVAL 
11.3    GAL/MAN/DAY    TO    20.0    GAL/MAN/DAY.     IT    GENERATES    DAILY 
OBSERVATIONS    OF    THE    DAILY    RATE    AND    CONVERTS    THEM    TO    HOURLY 
VARIATES    BASED    ON    AN    EMPIRICAL    ARRIVAL    DISTRIBUTION.     THE 
HOURLY    OBSERVATIONS    THEN    SERVE    AS    THE     INPUT    TO    A    SEWAGE 
HOLDING    TANK    FROM    WHICH    THE    SEWAGE    IS    PROCESSED    ACCORDING    TO 
ONE    OF    FIVE    SAMPLE    PROCESSING    POLICIES    WHICH    ARE    ARE    REPRE- 
SENTED   BY    INDIVIDUAL    SUBROUTINES.    OVER    A    WEEK'S    PERIOD    OF 
SIMULATED    TIME,     THE    MAXIMUM    SEWAGE    LEVEL    OBSERVED    FOR    UP    TO 
100    REPLICATES     IS    RECORDED    ALONG    WITH    THE    TOTAL    PROCESSING 
TIME.     STATISTICS    ARE    THEN    COMPUTED    AND    THE    THE    RESULTS    PRINT 
ED    ACCORDING    TO    ONE    OF    FOUR    POSSIBLE    PRINT    ROUTINES. 

JU  *V  -jU-  -ju  j.  «Js*  ***  *C  *V  V*  V*  **»  V"  *V  ***  V*  ***  */*  V*  V*  UV  *V  *V  V*  •*'*  ***  V*  ***  **•■  V*  §fr  "W  4f  Vr  *"*  *>**  **"  »**  »''  V^  **»  *V  ift*  *^*  i**  ***  ■***  "»V  **•  ^^  *A*  W*  ■UU  ■••*  >v  *V  -JU  si*  ***  JU> 

T  1-  i*  ',-  >,-  '.-  ,,.  T  ^  J,s  -,-  -r  -,-  t  t  i-  1'  1-  -i*  -t  t  'r  'i1  T  '»'  '.'  *,.-,-  *^  ^.  *r  T  -^  *,»  -r  ^,*  ^,.  -,.  *,,  ^,-.  -,,  ^  ^  ^  ^  -,.  ^,  ^  T  -^s  ^,  ^.  ^<  n,  ,,.  ^.  /,>  ,,v  r  r 

*  * 

THE    FOLLOWING    STATEMENTS     INITIALZE    THE    PROGRAM.    THE    DATA 
STATEMENT    READS    IN    THE    CONVERSION    FACTORS    FOR    THE    EMPIRICAL 
DISTRIBUTION 

COMMON    IX, I, J, ITIME,NOPRIN 

DIMENSION    CONVERT  24) ,X(7,24),Z(7?24) 

DIMENSION    DAILY (7) , HOURLY ( 7 ) , Y ( 7 ) ,  Yl ( 7) 

DIMENSION    X0RDQ68)  ,Y0RD(168) 

DIMENSION    ZMAX( 10, 10  0) , I  PROS S( 10,100) 

REAL*8    TITLE(12)/'TANKSIM2« , 11*'  •/ 

REAL    LABEL/'  •/ 

DATA    C0NVER/.756,1.97, 1.82,. 963,. 812, 1.09, 1.80, 1.61, .9 
D63,. 82 5, 1.47, 1.69, 1.31, 1.04,. 839,. 77  0, .880, .770,. 68 7, 
E.495, .344, .234, .440,.  413/ 

THE    PARAMETERS    TO    BE    TESTED    ARE    THEN    READ    IN 

RE AD ( 5, 20 ) MEN, NP LAN, ALPHA2 , ZSTOP ,Z STAR T, NOPR IN ,NODRAW 
CKSTOP 
20    F0RMAT(I3,I2,F4.1,F5.2,F5.1,3I2) 

THE    FIVE    POLICIES    ARE    THEN    TESTED    SEQUENTIALLY    FOR    100    REP- 
LICATIONS   EACH.    THE    SAME    RANDOM    NUMBER    SEED    IS    USED    FOR    EACH 
POLICY 

DO    9999    NRUN=1,4 
IX=21111111 
NPLAN=NPLAN+1 
ALPHA1=ALPHA2 

EACH    POLICY    IS    TESTED    FOR    NINE    RATES     IN    INCREMENTS    OF    10    GAL 

DO    6000    M=l,9 

ALPHA1=ALPHA1+10.0 

DO    5000    L=l,100 

IX=IX+112 

IF(N0PRIN.LT.2)     GO    TO    155 

WRITE(6, 150JL 
150    FORMAT (• 1' ,20X, 'RUN    NO.       «,I2,//) 
155    ZMAX{M,L)=0.0 

ALPHA=0.0 

ITIME=0 

DO    30    1=1,7 
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25 
30 

DO  25  J=l,24 
Z(If J)=0.0 

X( I tJ)  =  o.o 

CONTINUE 
CONTINUE 

SEVEN 
* 

BETA  VARIATES  ARE  GENERATED 

FOR  EACH  DAY  OF  THE  WEEK 

DO  50  1=1,7 
Y( I)=GAM(4.3,1.,1. ) 
YKI  )=GAM(5.7,1.,1.) 
BETA2=Y(I)/(Y(I  )+YKI)  ) 
DAILY ( I )=BETA2*8.7+11.3 
50   CONTINUE 

IF(N0PRIN.LT.2)  GO  TO  56 
WRITE (6, 55) (DAILY ( I), 1=1,7) 

55  FORMATCO  DAILY  BETA  VARIATES:  •  //  ( IX,  7  (  F5  .2  ,  IX  )  )  ) 

*  i 

THE  DAILY  VARIATES  ARE  CONVERTED  TO  HOURLY  VALUES  AND 
THEN  TRANSFORMED  TO  THE  EMPIRICAL  DISTRIBUTION 

$  3 

56  DO  60  1=1,5 

60  HOURLY(I)=(DAILY(I )*MEN)/24 

DO  65  1=6,7 
65  HOURLY(I)=(DAILY( I )*MEN/4)/24 

DO  60  1=1,7 

DO  70  J=l,24 

X(I , J)=CCNVER(J)*HOURLY( I) 
70  CONTINUE 
80  CONTINUE 

DO  1500  1=1,7 

IF(  (NPLAN.EQ.l) .OR. (NPLAN.EQ.4)  )  GO  TO  90 

*  j 

THIS  SECTION  THEN  TAKES  THE  SEWAGE  INPUT  THAT  HAS  BEEN  GEN- 
ERATED AND  PROCESSES  IT  ACCORDING  TO  ONE  OF  THE  FIVE  POLIC- 
IES. THE  RESULTS  FOR  EACH  HOUR  ARE  RECORDED  IN  A  MATRIX. 
THE  VALUE  FOR  THE  MAXIMUM  SEWAGE  LEVEL  IS  ALSO  RECORDED  FOR 
EACH  WEEKLY  CYCLE 

*  : 

ALPHA=ALPHA1 
90    DO    1000    J=l,24 

IF(  (  J.EQ.D.AND.  (I  .EQ.l)  )G0    TO    855 

IF(  (  J.EQ.D.AND.  (I  .GE.2)  )G0    TO    860 

Z(I,J)=    Z(  I,  J-D+XII  ,J)-ALPHA 

GO  TO  880 
855  IF  (NPLAN.EQ.4)  ALPHA=ALPHA1 

ZU,J)=X(I,J)-ALPHA 

GO  TO  880 
860  Z(I,J)=  ZU-1  ,24  )+X(  I,  JJ-ALPHA 
880  IF(Z( I ,J) .LT.0.0)Z(  I, J  )  =  0.0 

IF(ZMAX(M,L) .LT.Z(  I,  J)  )  ZMAX ( M, L )  =  Z( I , J ) 

TANK=Z(I,J) 

IF(NPLAN.EO.l)  GO  TO  1000 

IF(NPLAN.EQ.2)  GO  TO  1030 

IF(NPLAN.EQ.3)CALL  PLAN3 ( TANK, AL PHA1 , ALPHA , £900) 

IF (NPLAN.EQ.4) CALL  PLAN4(TANK, AL PHA1 , ALPHA, Z STOP , 
1ZSTART,&900) 

IF(NPLAN.EQ.5)  CALL  PL  AN5(  ALPHA  1 ,  ALPHA  ,  KSTOP  ,  £1900  ) 
900  CONTINUE 
1000  CONTINUE 
1500  CONTINUE 

:fe  : 

THE  PROCESSING  TIME  IS  ALSO  COMPUTED  AND  RECORDED 

&  2 

IF     (NPLAN.EQ.l)     ITIME=168 
IF(NPLAN.EQ.2)     ITIME=0 
IF(NPLAN.EQ.5)     IT  I ME= ( 24-KSTOP ) *7 
IPR0SS(M,L)=168-ITIME 
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J,  ,JL, 

IF  DESIRED  THE  WEEKLY  SEWAGE  LEVEL  CAN  BE  PLOTTED  FOR  EACH 
HOUR 

IF(NODRAW.EQ.O)  GO  TO  2900 

N=0 

DO  2000  1=1,7 

DO  1900  J=l,24 

N=N+1 

YORD(N)=Z( I, J) 
1900  CONTINUE 
2000  CONTINUE 

DO  2100  J=l,168 
2100  XORD(J)=J 

CALL  DRA W ( 168, XORD,YORD, 0,0, LABEL, TITLE, 2 0.0, 300. 0,0,0 
1LAST) 

THE  RESULTS  CAN  THEN  BE  PRINTED  TO  THE  EXTENT  DESIRED 

2900  IF(N0PRIN.LT.2)  GO  TO  2950 

WR IT E( 6,3500) ZMAX( M,L) , ALPHA 1,IPR0SS(M,L) 

3500  FORMAT  CO' ,20X,'ZMAX=' ,F8.1,10X, 'PROCESSING  RATE  =  », 
1F6.1,10X, "PROCESSING  TIME  =  ',13,//) 

2950  IF(N0PRIN.LT.3)  GO  TC  5000 

WRITE  (6,  3000)  UZ(I,J),I  =  1,7),J  =  1,24) 

3000  F0RMATC1   Z  MATRI  X'  //  (IX,  7  (  F8  .  1  ,4X  )  )  ) 
WRITE(6,3100)((X(I,J),I=1,7),J=1,24) 

3100  FORMATCO   X  MATRIX  ' // ( IX,  7  (  F8  .2  ,3X  )  )  ) 

5000  CONTINUE 

6000  CONTINUE 

THIS  ROUTINE  COMPUTES  THE  AVERAGE  VALUES  AND  VARIANCES  FOR 
THE  MAXIMUM  SEWAGE  LEVEL  AND  THE  WEEKLY  PROCESSING  TIME 

DIMENSION  ZSUM(20)  , ZA VE( 20 )  , ZVSUM( 20 )  ,  TSUM(20 ) , 
1TAVE(20) ,TVSUM(20) ,ZVAR(20) ,TVAR(20) 

DO  7200  J=1,M 

ZSUM( J)=0. 

ZVSUM( J)=0. 

TSUM( J)=0. 

TVSUM( J)=0. 

DO  7000  1=1, L 

ZSUM( J)=ZSUM(J)+ZMAX( J, I ) 

TSUM(J)=TSUM( J)+IPROSS(J,I ) 
7000  CONTINUE 

ZAVE( J)=ZSUM(J)/L 

TAVE(J )=TSUM(J)/L 

DO  7100  1=1, L 

ZVSUM( J)=ZVSUM(J)+( (ZMAX(J, I)-ZAVE( J) )**2) 

TV SUM ( J)  =  TVSUM(J )  + ( ( I  PRO SS( J , I )-TAVE( J) )**2) 

ZVAR( J)=ZVSUM( J)/L 

TVAR(J)=TVSUM( J)/L 
7100  CONTINUE 
7200  CONTINUE 

WRITE(6,8000)  NPLAN 
8000  FORMAT  CI     ZMAX  MATRIX  FOR  POLICY  NO.  ',12//) 

DO  8200  J=1,L 

WRITE(6,8100)(ZMAX(I ?J) ,I=1,M) 
8100  FORMAT (IX, 10(F8.2,3X) ) 
8200  CONTINUE 

WRITE (6, 8500) (ZAVE( J) ,J=1,M) 
8500  FORMATCO     AVERAGE  ZMAX'  //(  3X,  10(  F8.  2  ,  2X)  )  ) 

WRITE (6, 8600) ( ZV AR ( J )  ,  J  =  l ,  M  ) 
8600  FORMATCO     ZMAX  VARIANCE  '  //  (  3X,  1  0(  F  1  0.  2  ,  IX)  )  ) 

WRITE(6,9000) 
9000  FORMAT  CI    PROCESSING  TIME  MATRIX'//) 

DO  9200  J=1,L 

WRITE(6,8150)( IPROSS(I,J), I=1,M) 
8150  FORMAT (IX, 10(13, 3X) ) 
9200  CONTINUE 

WRITE (6,9300 )(TAVE( J ) ,J=1,M) 
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9300  F0RMATP0     AVERAGE  TIME*  //(3X»  10  (  F6.2,2X)  ) ) 

WRITE (6, 9400 )(TVAR( J  ) , J=1,M) 
9400  FORMAT  CO     TIME  VARIANCE  • // ( 3X, 10 ( F 10. 2 , IX) ) ) 
9999  CONTINUE 

STOP 

END 

FUNCTION  GAMULFA, BETA, START) 
C  NON  INTEGER  GAMMA  GENERATOR  DEVELOPED  BY  D  T  PHILLIPS 
COMMON  IXtltJf ITIMEi NOPRIN 
IF( START. GT.l. 5)  GO  TO  58 
X3=1.0 

IFCALFA.LE.2.0)  GO  TO  1 
IFULFA.LE.5.0)  GO  TO  2 
GO  TO  3 

1  B=0. 2 177+ (2. 10299768* ALFA) -(4. 3461 196 1*ALFA**2)+ 
C(6c00914764*ALFA** 
C3)-(3.95097728*ALFA**4)+(0.97279251*ALFA**5) 

GO  TO  4 

2  B= 0.643 50+ (0.45839602* ALFA )-(0. 02952 801*ALFA**2)+ 
C(0.00172718*ALFA* 

C*3)-( 0.00005 8 10* ALFA**4)+( 0. 000000 8 2* ALFA** 5 ) 
GO  TO  4 

3  B=l. 33 408 +( 0. 2249999 1*ALF A )-(0. 0023069 5*ALFA**2)+ 
C( 0.00001623*ALFA* 

C*3)-(0. 00000006* ALFA**4) 

4  S=1.0+(1.0/B) 
START=5.0 

12  IF(S-1. 0)110,60, 15 
15  S=S-1.0 

X3=X3*S 

GO  TO  12 
110  GY=1.0+S*(-0.5771017+S*(0.98  5  854+S*(-0.8  764218+S*(0.83 
lS*(-0.56  84729+S*(0.2  54  8205+S*(-0.0  5149930) )))))) 

X3=X3*GY 
58  A=(X3/(ALFA*BETA) )**B 

B=1.0/B 

A=1.0/A 
60  CALL  RANDU(IX,IY,YFL) 

IX=IY 

GAM=(-A*ALOG(YFL) )**B 

RETURN 

END 

SUBROUTINE  PLAN3 (T ANK, ALFA1 ,RATE,*) 

THIS  POLICY  PROCESSES  SEWAGE  UNTIL  1600  AFTER  WHICH  TIME  IT 
THEN  EXAMINES  THE  LEVEL  OF  THE  TANK  AND  IF  BELOW  50  GAL. IT 
STOPS  PROCESSING  UNTIL  0500  THE  FOLLOWING  DAY.  ON  WEEKENDS 
PROCESSING  STOPS  ANYTIME  THE  LEVEL  IN  THE  TANK  REACHES  50  GA 

COMMON  IX, I, J, ITIME, NOPRIN 

IF( ( I .EQ.l) .AND. (J.EQ.l) )  IFLAG=0 

IF(N0PRIN.LT.2)  GO  TO  6 

IF( ( I . EQ.l ).AND.( J.EQ.l) )  WRITE (6, 5) 

5  FORMAT  CO     PROCESSING  POLICY   NO. 3  WAS  USED  ON  THIS 
CRUN' ,//) 

6  IF( (J.EQ.l) .AND. (IFLAG. EQ.l) )  I T IM E= ITI ME  +  24- JSTOP 
IF(J.EQ.l)  IFLAG=0 

IFU.LE.5)  GO  TO  10 

IFUANK.LE.50.  )  GO  TO  30 

IF( (RATE. LT. 0.01). AND. (TANK. GE. 500.) )  GO  TO  20 

GO  TO  40 
10  IF(J.LE.ll)  GO  TO  40 

IF(TANK.LE.50. )  GO  TO  30 

GO  TO  40 
20  RATE=ALFA1 

IF(IFLAG.EQ.O)  GO  TO  40 

ITIME=ITIME+J-JSTOP 

IF(N0PRIN.LT.2)  GO  TO  26 
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WRITE(6,25)I ,J 

25  FORMAT P     CONTINUED  PROCESSING  AT  T= • , 1 1  ,  •- '  ,  I  2 , // ) 

26  IFLAG=0 
GO  TO  40 

30  RATE=0.0 

IFCIFLAG.EQ.DGO    TO    40 
JSTOP=J 

IF(N0PRIN.LT.2)     GO    TO    36 
WRITE(6,35)I ,J 

35  FORMATt'O      STOPPED  PROCESSING  AT  T= • , 1 1 , •-■ , I  2 , // ) 

36  IFLAG=1 

40  IF( ( I .EQ.7).AND.( J.EQ.24) )  GO  TO  45 

GO  TO  50 
45  IF(IFLAG.EQ.l)  IT IME=IT IME+24-JST0P 
50  RETURN  1 

END 

SUBROUTINE  PL AN4 (TANK, ALFA  1 , RATE , STOP, ST ART f* ) 

THIS  IS  A  LEVEL  CONTROL  POLICY  WHICH  CONTINUOUSLY  EXAMINES 
THE  LEVEL  OF  SEWAGE  IN  THE  TANK. IF  BELOW  SOME  SET  LEVEL 
(ZSTOP)  THEN  PROCESSING  STOPS  UNTIL  THE  SEWAGE  LEVEL  REACHES 
SOME  UPPER  VALUE  (ZSTART)  AT  WHICH  TIME  PROCESSING  CONTINUES 

COMMON  IX,  I, J, ITIME,NOPRIN 

IF(  (  I.EQ.D.AND.  (J.EQ.l)  )  IFLAG=0 

IFCNOPRIN.LT. 2)  GO  TO  6 

IF((  I.EQ.D.AND.  (J.EQ.l)  )  WRITE  (6,  5) 

5  FORMAT («0     PROCESSING  POLICY   NO. 4  WAS  USED  ON  THIS 
CRUN' ,//) 

6  IF( (J.EQ.l). AND. (IFLAG.EQ. 1) )  I T IM E= ITIME+24-JST0P 
IF(J.EQ.l)  JSTOP=0 

IF(TANK.GT.STOP)  GO  TO  20 

RATE=0.0 

IFdFLAG.EQ.DGO    TO    40 

JSTOP=J 

IF(N0PRIN.LT.2)  GO  TO  11 

WRITE( 6,10) I , J 

10  FORMATt'O      STOPPED  PROCESSING  AT  T=" , 1 1  ,  »-• , 12, //) 

11  IFLAG=1 
GO  TO  40 

20  IF(TANK.LE. START)  GO  TO  40 

RATE=ALFA1 

IF( IFLAG.EQ. 0)  GO  TO  40 

ITIME=ITIME+J-JSTOP 

IF(N0PRIN.LT.2)  GO  TO  35 

WRITE(6,30)I ,J 
30  FORMATC     CONTINUED  PROCESSING  AT  T=  •  ,  II  ,  ■-•  ,  I  2,  //) 
35  IFLAG=0 
40  IF( ( I.EQ.7).AND. (J.EQ.24) )  GO  TO  45 

GO  TO  50 
45  IF(IFLAG.EQ.l)  I T I ME=ITIME+24-JST0P 
50  RETURN1 

END 

SUBROUTINE  PLAN5 ( ALF Al , RATE , KTI ME, *) 
&  $ 

THIS  IS  A  TIME  CONTROL  POLICY.  IT  ALLOWS  PROCESSING  FOR  A 

FIXED  PORTION  OF  THE  DAY  ONLY 

$  $ 

COMMON  IX,I,J,ITIME,NOPRIN 

IF(N0PRIN.LT.2)  GO  TO  6 

IF(  (  I.EQ.D.AND.  (J.EQ.l)  )  WRITE  (6,  5) 

5  FORMATCO     PROCESSING  POLICY   NO. 5  WAS  USED  ON  THIS 
CRUN' ,//) 

6  IF(J.LE.KTIME)  GO  TO  10 
RATE=0.0 

10  RETURN  1 
END 
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indicated  that  a  revision  of  the  Navy's  design  parameter  for  the 
dailv  per  capita  sewage  generation  rate  was  in  order.   The  Simula 
Uons  were  aesigned  so  that  use  by  others  with  their  own  data  would 
be  an  easy  matter. 
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